HOPF FIBRATION: GEODESICS AND DISTANCES 
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Abstract. Here we study geodesies connecting two given points on odd-dimensional 
spheres respecting the Hopf fibration. This geodesic boundary value problem is completely 
solved in the case of 3-dimensional sphere and some partial results are obtained in the 
general case. The Carnot-Caratheodory distnee is calculated. We also present some 
motivations related to quantum mechanics. 



1. Introduction 

Sub-Riemannian geometry is proved to play an important role in many applications, e.g., 
in mathematical physics, geometric mechanics, robotics, tomography, neurosystems, and 
control theory. Sub-Riemannian geometry enjoys major differences from the Riemannian 
being a generalization of the latter at the same time, e.g., the notion of geodesic and length 
minimizer do not coincide even locally, the Hausdorff dimension is larger than the manifold 
topological dimension, the exponential map is never a local diffeomorphism. There exists a 
large amount of literature developing sub- Riemannian geometry. Typical general references 
are P H El 120]. 

The interest to odd- dimensional spheres comes first of all from finite dimensional quan- 
tum mechanics modeled over the Hilbert space C n where the dimension n is the number 
of energy levels and the normalized state vectors form the sphere § 2 ™ -1 C C n . The prob- 
lem of controlled quantum systems is basically the problem of controlled spin systems, 
which is reduced to the left- or right-invariant control problem on the Lie group SU(n). 
In other words, these are problems of describing the sub- Riemannian structure of S 2n_1 
and the sub- Riemannian geodesies, see e.g., jH [15]. The special case n = 2 is well stud- 
ied and the sub- Riemannian structure is related to the classical Hopf fibration, see, e.g., 
[2"Ij |2"2"] . At the same time, the sub- Riemannian structure of § 3 comes naturally from the 
non-commutative group structure of SU(2) in the sense that two vector fields span the 
smoothly varying distribution of the tangent bundle, and their commutator generates the 
missing direction. The missing direction coincides with the Hopf vector field corresponding 
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to the Hopf fibration. The sub-Riemannian geometry on § 3 was studied in [HI [TTJ [13] , see 
also (5J. Explicit formulas for geodesies were obtained in [9J by solving the corresponding 
Hamiltonian system, in [TT] from a variational equation, in [5] exploiting the Lie theory, in 
[H] applying the structure of the principle § 1 -bundle. One of the important helping prop- 
erties of odd- dimensional spheres is that there always exists at least one globally defined 
non- vanishing vector field. 

Observe that S 3 is compact and many properties and results of sub-Riemannian geometry 
differ from the standard nilpotent case, e.g., Heisenberg group or Engel group. In the case 
§ 2n_1 , n > 2, we have no group structure and the main tool is the global action of the 
group U(l). For example, in our paper we explicitly show that any two points of S 3 can 
be connected with an infinite number of geodesies. 

Because of important applications, we start our paper with the description of n-level 
quantum systems and motivation given by Berry phases. Further we continue with general 
formulas for geodesies. Then we concentrate our attention on the geodesic boundary value 
problem finding all sub-Riemannian geodesies between two given points. In the case of 
§ 2n_1 we solve it for the points of the fiber and for S 3 we solve it for arbitrary two points. 
The Carnot-Caratheodory distance is calculated. 

Acknowledgement. This work was initiated while the authors visited the National Cen- 
ter for Theoretical Sciences, Hsinchu, Taiwan R.O.C. The authors acknowledge support and 
hospitality extended to them by the director of the center professor Wen-Ching Winnie Li 
and the staff. 



2. n-LEVEL QUANTUM SYSTEMS 

The mathematical formulation of quantum mechanics is based on concepts of pure and 
mixed states. A complex separable Hilbert space $) with Hermitian product (•, •) (or (-|-) in 
Dirac notations) is called the state space. The exact nature of this Hilbert space depends on 
the concrete system. For an n-level quantum system, Sj = C n with the standard Hermitian 
product (z\w) = Yjj=i z j™j- An observable is a self-adjoint linear operator acting on the 
state space. A state p is a special case of observable which is Hermitian pt = p, normalized 
by tr p = 1, and positive (w\p\w) > for all vector \w) G S), where p\w) denotes the result of 
the action of the operator p on the vector \w). A pure state is the one-complex-dimensional 
projection operator p = \z)(z\ onto the vector \z) G $], i.e., an operator satisfying p 2 = p. 
Other states are called mixed states. 

The space of pure states is isomorphic to the projectivization i^p of the Hilbert space 
S). So equivalently we can define pure states as normalized vectors (z\z) = 1 modulo a 
complex scalar e te , where 9 is called a phase. In the case of the n-level quantum system 
Sj = C n , normalization and the phase factor allow us to represent the space of pure states 
as the complex projective space CP" -1 = Cp. The second operation of phase factorization 
is realized by the higher-dimensional Hopf fibration 

gl ^ gJJn-1 _^ C¥ n-1 
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where S 2n_1 is the real (2n — l)-dimensional sphere embedded into C n , and the base space 
CP™ -1 is the set of orbits of the action S 1 on S 2n_1 . 

In what follows, the pure states p are the elements of the projective space CP n_1 , at the 
same time we use the notation of state \z) G S 2 ™ -1 to refer to a normalized representative 
of the phase-equivalence class p. The real dimension of the space of pure states is 2n — 2. 
The projective space CP™ -1 endowed with the Kahlerian Fubini-Study metric on CP n_1 
becomes a metric space, in which this Riemannian metric is given by the real part of the 
Fubini-Study metric and it coincides with the push-forward of the standard Riemannian 
metric (the real part of the Hermitian one) on the (2n — l)-sphere by the Hopf projection. 

In the case n = 2 the projective complex plane CP is isomorphic to the sphere S 2 (called 
the Bloch sphere for 2-level systems in physics), which is thought of as the set of orbits of 
the classical Hopf fibration 

S 1 § 3 ->■ S 2 . 

Each pair of antipodal points on the Bloch sphere corresponds to a mutually exclusive pair 
of states of the particle, namely, spin up and spin down. The Bloch sphere and the Hopf 
fibration describe the topological structure of a quantum mechanical two-level system, see 
[2T| [22]. The interest to two-level systems, an old subject, recently has gained a renewed 
interest due to recent progress in quantum information theory and quantum computation, 
where two-level quantum systems became qubits coupled in g-registers. A qubit state is 
represented up to its phase by a point on the Bloch sphere. The topology of a pair of 
entangled two-level systems is given by the Hopf fibration 

§ 3 ^ § 7 HP, 

where H is the space of quaternions and HP = § 4 is its projectivization, see [TH]. Generally, 
for entangled n- level systems we have 

The underlying manifold for the Lie group SU(2) is S 3 . Considering the higher dimensional 
group SU(n), we see that it acts on § 2n_1 . However its dimension is n 2 — 1 > 2n — 1, 
n > 2, and its manifold only contains the invariant sphere S 2n_1 . Returning back to the 
information theory motivation, the relevant group for p-qubits is SU(2 P ), see, e.g., [6| [T8]. 

Let us now concentrate on n-level systems. Time evolution of a quantum system in the 
Schrodinger picture assumes time-dependent states while observables are conserved. Time 
evolution of a quantum state is described by the Schrodinger equation 

j t \m = ^H\ Z (t)), 

where if is a special observable, the Hamiltonian, a Hermitian operator corresponding to 
the total energy of the system. Its spectrum is the set of possible outcomes of energy 
measurements which we suppose to be discrete and non-degenerate. The normalization 
condition is preserved by time evolution of the quantum system. The Bloch equation for 
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the projection operator corresponding to the state \z(t)) becomes the Hamiltonian equation 



where [H,p] stands for the operator commutator. Observe that the state \z) is viewed 
as an element of § 2n-1 while the operator p is thought of as an element of CP™" 1 . Thus 
the solutions to the Schrodinger equation describe trajectories on § 2n_1 , and the solution 
to the Bloch equation describe trajectories on CP" -1 . The solution to the Schrodinger 
equation with the Hamiltonian H(t), —iH(t) G su(n), is of the form \z{t)) = g(t)\z(0)), 
where g(t) G SU(ra). Observe that if H is an Hermitian operator and if we denote by 
h = i(tri7 — 1), then H = H — h l is an Hermitian operator and tv(H — h l) = 1. 
The operator H = h Q l is called the free Hamiltonian. If the system is considered under 
some external fields (e.g., an electromagnetic field) whose amplitudes are represented by 
some functions hi(t), h m (t), then H(t) = H + Y^h=i hk(t)H k , where H fc are self-adjoint 
operators representing the coupling between the system and the external fields. One of 
the topical problems is population transfer, to find the amplitudes hk(t), k = 1 . . . m, such 
that the system starts at an eigenstate \zq) of the free Hamiltonian Hq and ends at time 
T at another eigenstate \z{) of H 

The controllability problem resides in the so-called bracket-generating condition by Chow 
and Rashevsky JTUJ [19] . The system is controllable if and only if the Lie hull 



Suppose that the operator pit), being viewed as an element of the projective space, 
evolves as a loop in time t G [0, T], p(0) = p(T). A crucial concept in many quantum 
mechanical effects is so-called the Berry phase (or Berry-Pancharatnam following [31 HT]). 
The state \z(t)) of the system evolves according to the Schrodinger equation. Suppose that 
\z{T)) = e i0 |z(O)). Then 6 = 6 geom + Q dyn , where O 9eom is the Berry phase and Q dyn 
stands for the standard dynamical phase given by 



Observe that the instantaneous energy (z(t)\H\z(t)) is real because H is self-adjoint. The 
geometric Berry phase refers to the curvature of the bundle. It emphasizes that the single- 
time probabilistic description does not exhaust the physical content of quantum mechanics 
and appears because of the bundle structure of the space of states (see, e.g., [2]). As 
Berry [3] mentions, Barry Simon commented in 1983 that the geometric phase factor has 
an interpretation in terms of holonomy because the phase curvature two-form emerges 
naturally as the first Chern class of a Hermitian line bundle. 

Let us suppose that the Hamiltonian H is perturbed H £ u\ by a varying parameter e(t), 
e(0) = s{T) slowly enough so that the adiabatic approximation is applied with T large. 
The systems undergoes transport along the closed path in the parametric space e. The 
adiabatic theorem states that if the system starts at an eigenstate of the Hamiltonian 
if e (o), then it will end in the corresponding eigenstate of the Hamiltonian H e (t) at the final 
point under small perturbation remaining in its instantaneous eigenstate in all intermediate 



Lie (iH , iHi, . . . , iH m ) = su(n). 
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points t G [0, T\. This theorem has been known from early years of quantum mechanics 
due to Max Born and Vladimir Fock. 

Once the controllability problem is solved one can look for a more efficient way to 
construct loops in the projective space of pure states. A typical way is to minimize the 
energy of transfer. The corresponding optimal control problem is equivalent to finding 
sub-Riemannian geodesies on S 2n_1 . The precise meaning of this will be explained in the 
following section. 

3. SUB-RlEMANNIAN STRUCTURE AND GEODESICS ON § 2 ™ -1 

3.1. Hopf fibration. We consider the (2n — l)-dimensional sphere § 2n_1 as a real manifold 
embedded into the complex space C n . The Riemannian metric d on the sphere is the 
restriction of the Euclidean metric from R 2n , that is, the real part of the Hermitian product 
(•, •) in C n . The complex Hopf fibration h is a principal circle bundle 



over the complex projective space CP n_1 . This roughly means that the preimage h~ l (m) 
of any point m G CP" -1 is isomorphic to the unit circle S 1 endowed with the Lie group 
structure of U(l). The set h~ l {m) C § 2n_1 is called the fiber, and CP 71-1 is the base space, 
see, e.g., 0Q3]. 

The intersection of a complex line with the unit sphere S 2n_1 gives a unit circle S 1 that 
is equipped with the group structure of U(l). The action of e 2mt G U(l) on § 2n_1 is defined 
by 



We will use both real (xi, yi, . . . , x n , y n ) and complex coordinates z k = x k + iyk, k = 
l,...,n. 

The horizontal tangent space T-LT Z at z G S 2n_1 is the maximal complex subspace of the 
real tangent space T z . The unit normal real vector field N(z) at z G S 2 " -1 is given by 



k=l k=l 



N(z) = ^2(x k d Xk + y k d yk ) = 2Re ^ 

The real vector field 

V(z) = iN(z) = ^2(-y k d Xk + x k d yk ) = 2Re ^ iz k d. 



k=l k=l 



is globally defined and non-vanishing and we call it the vertical vector field. The vector 
space span M {iV(,2), V(z)} is a two dimensional subspace of C n that inherits the standard 
almost complex structure from C n , in other words we have an isomorphism 

spanJiV^), V(z)} = span c {iV(z) + iN(z)} 

of vector spaces. The orthogonal complement to span c {iV(z) + iN(z)} with respect to 
the Hermitian product is exactly the horizontal tangent space "HT 2 C T z . If we denote 
by ip(t) = e 27Tlt z, t G [0,1), the fiber at a point z G S 2n_1 , then it is easy to see that 



6 



DER-CHEN CHANG, IRINA MARKINA, AND ALEXANDER VASIL'EV 



(p{t) = 2irV(<p(t)). Thus, span M {V(£)} is the tangent space of the fiber at z = <p(t), that 
we denote by TT Z . The vertical tangent bundle of <p we denote by TT, and the horizontal 
tangent bundle is denoted by %T. Then, 

(1) T Z = HT Z ®TT Z , 

where the direct sum means the orthogonal decomposition with respect to the Riemannian 
metric d at each z G 8 2n_1 . The restriction of d to the horizontal and vertical subspaces is 
denoted by d n and djr, respectively. 

The group SU(n) acting in C n leaves S 2n_1 invariant and preserves the orthogonal de- 
composition (frj. 

Concluding the above construction we will work with a sub-Riemannian manifold that 
is the triple (S 2n_1 , TiT, dy) consisting of the smooth manifold § 2n_1 , the horizontal sub- 
bundle (or distribution HT), and the sub-Riemannian metric dy_. The same horizontal 
bundle can be obtained by considering the standard C-R geometry of odd-dimensional 
spheres or thinking of § 2n_1 as a contact manifold endowed with the one- form [131 EH] 

n 

(2) n = y^X-Uk dxk + %k dVk)- 

k=i 

We are looking for sub-Riemannian geodesies, or energy minimizing curves j(t) on S 2n_1 . 
Since the curves have to belong to the sphere and to be horizontal, we conclude that the 
tangent vector j(t) is orthogonal to complex line span c (iV + iN). Indeed, 

(3) Re ( 7 , 7> = 0, Re (z 7 , 7) = "I™ <7, 7> = (7, 7) = 0, V t € [0, 1]. 



3.2. Equations of geodesies. The Hamiltonian function in sub-Riemanian geometry is 
defined by using the notion of co- metric [151 120] • As it was mentioned above, the sub- 
Riemannian metric dy_ is defined by the restriction of the Euclidean metric from R 2n to 
the horizontal sub-bundle HT. The restriction of the Euclidean metric from R 2n to fibers 
TT gives the metric djr with respect of which the vertical vector field V has length one. 
Given d^ we may define a linear mapping g%(z) : T* — » T z . This linear map is uniquely 
defined by the following conditions: 

1. im.(g H (z)) = 1-LT Z at any point z G § 271_1 , 

2. p(v) = d H (g H (z)\p\,v) for v G UT Z , p G T* at each z G S 2 ™" 1 . 

We denote by Tj~ the one-dimensional kernel of gu{ z ) at 2 G S 2 ™ -1 . The generator of the 
kernel is the contact form or connection one- form (T5]). The linear map gu{z) is positively 
definite on T* z jT z and symmetric: p(g H (z)[u]) = d n (g n (z)\p), g n {z)[u\) = u(g H (z)\p\), 
p,u G T*, z G § 2n_1 . The map g% is called co-metric. 
Analogously, we define a linear map gr(z) : T* — >■ T z by 

1. im(gjr(z)) = TT Z at any point z G § 2n_1 

2. p(w) = d F {g T {z)\p\,w) for w G TT Z , p G T* at each z G S 2 " -1 . 
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Now we are able to set up two Hamiltonian functions: the horizontal Hamiltonian func- 
tion H n (z,p) and the vertical one Hj:(z,p), (z,p) G T*. We write 

Hu(z,p) = Ip(9h(z)\p\), H T (z,p) = ]zp{9f{z)\p\). 

Lemma 1. The horizontal H n and vertical Hjr Hamiltonian functions are Poisson invo- 
lutive: 

{H n ,H T } = 0. 

Proof. Consider a neighborhood in T*^ 2 ™" 1 constructed by means of a local trivialization 
of the circle bundle 

Let U be an open neighborhood in CP n_1 , W = h~ l (U), and ip is a local trivialization 
map ip : W — >■ U x S 1 . The linear map ip*{z) generated by ijj defines a diffeomorphism on 
cotangent bundles T*W ^T*U x T*§ x . 

Notice that the decomposition TU x TS 1 is orthogonal with respect to the metric d, and 
the metrics d^ and rfj- were defined as the restriction of d onto TU and TS 1 respectively. 
We write = (x 1 , . . . ,x 2n ~ 2 ,pi, . . . ,P2n-2) £ T*[7 and (y, /x) G T*S X within the proof 

of the theorem. Thus, the metric djr has unit real 1 x 1-matrix in the basis of the vertical 
vector fields V and the corresponding co-metric gjr defined by the same matrix. Therefore, 
they are independent of the points z G S 2ri_1 . Moreover, since V is normal and orthogonal 
to T-L everywhere we conclude that Hj?(x,p,y, p) = \p? 1 where p is the dual form to V 
and counting the dimension of p G T- 1 . The metric d-^ is independent of y by definition 
so does the co- metric gu- Moreover, the co- metric gu is independent of p by p G T- 1 . We 
conclude, that H n depends only on (x,p). We calculate 

2 "~ 2 'dH H d¥L T dH H dH F \ dE H dE T dH n dH T 



1 ^ V dvb dx k dx k dvb J 



0. 



V dpk dx k dx k dpk ) dp dy dy dp 

□ 

Theorem 1. A normal sub-Riemanian geodesic •y(s) starting at a point a G § 2n_1 is given 
by the formula 

(4) 7 (s) = lR (s)e- isd(v ' v{a)) = (acos(\\v\\s) + sm(\\v\\s))e- isd{v ' via)) , 

\\v\\ 

where v is the initial velocity of the Riemannian geodesic Jr(s) : 7fl(a) = v at the initial 
point a, and d(v, V(a)) is the projection of the initial velocity to the vertical direction at a. 

Proof. We start with a discussion of relations between two definitions of Riemannian 
geodesies: by the Riemannian exponential map and as the projection of the solution to 
the Hamiltonian system to the underlying manifold. It is a well known connection, but 
we recall it for completeness. To obtain a Riemannian geodesic Jr(s) starting at a point 
a G S 2 ^ 1 with the initial velocity vector v G T a S> 2n_1 we can make use of the Riemann- 
ian exponential map exp d (sv ) defined by the metric d. Let us recall, that given a metric 
d = {dtj} the Riemannian Hamiltonian function H d (z,p) = \d % ipi{z)pj{z) is defned, where 
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the co-metric d^ is the inverse matrix for the metric d. The solution to the Hamilton- 
ian system defines the flow &(q ,p ,t) on the cotangent bundle T*S 2n_1 , where (qo,Po) is 
the initial point. Moreover, the metric tensor d defines the identification of tangent and 
cotangent bundles. Thus, the Riemannian geodesic 'Jr(s) = exp d (sv) is the result of the 
following superpositions 

exp d (sv(a)) : T a § 2 " -1 A TS 2 ^ 1 ^ T*^ 1 A S 2 ™" 1 . 

Here the first map i is the inclusion of the tangent space at a G S 2n_1 into the tangent 
bundle T*§ 2n_1 . The second map is the identification of the tangent and cotangent bundles 
by means of the metric d, which also gives the initial data (qo,po) = (o, d(v, •)) for the flow 
$(go,Po, s). The last map projects the Hamiltonian flow $(go,Po, s) to the sphere. 

Pick up a vector v G TJs 2n ~ l and write v = v-^ + vjr due to the orthogonal decom- 
position (JT[) . We can construct exp d (si^) and exp d ^(st>jr) by means of the procedure 
described above using the horizontal Hy_ and the vertical Hjr Hamiltonian functions. The 
identification between the tangent and cotangent spaces is realized by the linear maps 
d-n : HT —> T*S 2n ~\ g H : T*§ 2n_1 ->■ TS 2 "" 1 and d T : TT -> T*$ 2n -\ g T : T*^ n ~ x -> 
T§ 2n_1 , respectively. Then 

exp d {tv) = exp dn {tv n ) o exp^(t^) 

since the Hamiltonian functions H-n, Hjr, Hd are all Poisson involutive, the corresponding 
flows commute. We observe also that the group exponential map exp § i : u(l) — >■ U{1) 
agrees with the Riemannian exponential map exp djr : u(l) = TT a — > S 1 [12J. If we write 
j R (s) = exp d (sv), ~/ sR (s) = exp dn (sv H ), e is ^ = exp djr (svjr), then 

(5) j sR (t) = Ms)e~ lsM = Ms)e- tsd(v ' V{a) \ 

here we used the notation |t>jr| = d(v,V(a)) for the projection of the initial velocity onto 
the vertical direction V(a) at a. 

It is well known that the Riemannian geodesic on sphere is a great circle and can be 
written as 

v 

(6) 7r( s ) — acos(||u||s) + |j— jj-sin(||t> \\s). 

Here a is the starting point, v is the initial velocity, and ||t>|| 2 = d(v,v) is the real norm 
of v. Combining fl3j) and © we obtain (|3|. 

To be sure that the geodesic is horizontal we verify the product 

(isRilsR) = V(o))(7fl,7fl) + (7b, 7b). 

Since 7^ G S 2n_1 , then (7^, 7^) = 1, Re (tr, 7^) = 0, and we need to calculate Im ( A f R , 7^). 
We write a k = a k + i(5 k and v k = v k + iu k . Then V(a) = Y2=i(~f^kd Xk + oc k d yk ) and 

zlm (j R {s), 7b(s)} = «Im XlLi ( _ + ^)IMI sin (IMI s ) + ( v h + i^k) cos(||v||s) J 
x ((a, - ift) cos(||v||s) + sin(||«||s)) 
= *Z)fc=i {- PkVk + oikUk) =id(v,V(a)). 



HOPF FIBRATION. 



9 



This completes the proof of the theorem. □ 

Remark 1. We have used the ideas developed in [15J for general principal G-bundles, in 
our case for the circle bundle of (2n — 1)- dimensional spheres and the sub-Riemannian 
metric given by the restriction of the Euclidean metric. The equations for geodesies in 
the particular case of § 3 were obtained in [9] by solving the Hamiltonian system, in [5] by 
applying the general Lie group theory and in [11] as a solution to a variational equation. It 
was shown in [13] that sub-Riemannian geometry on § 3 induced by CR-stucture ofE> 3 ^ C 2 , 
by the Hopf fibration, and by Lie group structure § 3 of unit quaternions coincide. 

Remark 2. The solution to a variational equation in [UJ depends on a real parameter 
A that the authors call curvature. This parameter for a geodesic starting at the point 
a coincides with d(v,V(a)) as was shown in [14]. The variational equation can not be 
generalized into the higher dimensions, but as we will see the value d(v, V(a)) will play an 
important role in the behavior of geodesies. We continue to use the notation \vj?(a)\ = 
d(v,V(a)). 

Remark 3. The group SU(ra) acts in C™ transitively on § 2n_1 = SU(n)/ SU(n — 1), more- 
over, the action is conformal isometrically. Therefore, the normal vector field N is in- 
variant under the action of SU(n). Since the orthogonal structure %T © span c (iV + iN) 
is preserved under the action of SU(ra) we conclude that the vector field V = iN is also 
invariant under the action o/SU(n). Thus if £ SU(n), then 

TT a 3 V{a) H- faV(a) = V{<j)a) £ TT^, UT a ^ UT^. 

Furthermore, 

\vjr(a)\ = d(v, V(a)) = d((f>v, V(<t>a)) = \vj?((j)a)\, d(v, v) = d((f>v, <fw) 

v H (a) ^ <frv H (a) = v H ((f)a). 

Any geodesic parametrized by arc length, that is \\vu\\ = 1, is uniquely determined by 
the initial point a, and the initial horizontal v-n(a) and the vertical vjr(a) velocities. Let us 
denote such a curve by 7(5; a,v-u,vjr) . Then 

070; a, v-h, vj) = 7(5; 0a, ^v H , <j)*v T ). 

4. Open and closed sub-Riemannian geodesics on S 2 "^ 1 

Let us suppose that geodesic (TjJ is parametrized by arc length: = 1. Then ||t>|| = 

yjl + \v jr| 2 . We are looking for the points of intersection with the vertical fiber at a, whose 
equation is given by ae^, (ft £ [0, 2%). 

Theorem 2. Under the above notations we claim 

1. If j=d = £ < I is rational with p, q > relatively prime, then for n = (2q)m, 

m = 1,2, . . ., the geodesic is a closed curve which meets the base point a for the first time for 
m = l, and then, it is periodic with minimal period 2t\\/ q 2 — p 2 . The length of each closed 



10 DER-CHEN CHANG, IRINA MARKINA, AND ALEXANDER VASIL'EV 

/ * 

loop of the geodesic is 2it\J q 2 — p 2 . The geodesic intersects the fiber at points ae 
n = 1,2, ... , 2q, for the first time and then periodically returns back. The length of the 

I 2 2 

part of the geodesic between the intersection points is equal to 7Ti ' 9 p 



2. If , = is irrational, then the geodesic is open and is diffeomorphic to the straight 



inn 1 



\V£\ 



line. It meets the fiber for s = , = n, n = 1, 2, . . at the points ae v /l+l "- F ' z 

Proof. In order to proof the theorem we need to find the moment s — s, such that the 
equation 

e~ isM (a cos (sy/l + M 2 ) + ^ = sin ( Sy /l + \v r \ 2 ) ) = a 

is verified. Due to Remark [3] it is sufficient to check the base point a = (1, 0, . . . , 0). It 
follows that for s = s = , n n, n = 1,2, . . ., the geodesic intersects the fiber at the 



inn 1 



points ae V V 1+ i"^ 2 



I".fI 



«7rn 1 



If the geodesic is closed, then ae V 1+l "^ 12 = a. We conclude that if , 

is rational and -, = = 2 < 1 where p, g are relatively prime and positive, then for 

n = 2q the geodesic returns back to a for the fist time, and then it becomes periodic 
passing through the point a for each n = (2q)m, m = 1,2,3, . . .. The minimal period is 
obtained by setting n = 2q and \v j-| 2 = q 2_ p 2 m s and equals 2tx\/ q 2 — p 2 . The Carnot- 
Caratheodory length of one closed loop is equal to s = 2tx^/ q 2 — p 2 . The first loop of 
a closed geodesic at the moments s = — n, n = 1, 2, . . . , 2q, intersects the vertical 



circle at the points ae v 9 ' . The Carnot-Caratheodory length of each part between 



intersections is 



If the fraction -. = is irrational, then the geodesic is open and is diffeomorphic 

yi+kH 2 

to the straight line. It also intersects the vertical circle at the moments s = , = n 



inn 1 



at the points ae a/i+i^i 2 / . The Carnot-Caratheodory length of each part between 
intersections is , w □ 
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We can prove the following theorem (see also [11]) using the group structure of § 3 , which 
is isomorphic to SU(2). The elements of SU(2) are the matrices 



z\ z 2 
-z 2 Z\ 



1 



where zi, z 2 G C, and \zi\ 2 + \z 2 \ 2 = 1, hence we also denote the elements of SU(2) by unit 
vectors (zi, z 2 ). Denote by e the unity (1, 0) of the group. 

Theorem 3. The left translation 07 of a geodesic 7(s; e, Vy_, \vj:\) by G SU(2) ; = 
(p, iv y/l — p 2 ), where p 2 = (l + (y/l + \vjr\ 2 — |fjr|) 2 ) , belongs to the Clifford torus 
{(W!,w 2 ) G S 3 : \ Wl \ 2 = p 2 }. 

Proof. To prove the theorem we need to find an element = (0i, <p 2 ) G SU(2), such that 

(pj(s;e,v H , \vjr\) = (wi(s),w 2 (s)), with \ Wl \ 2 = p 2 , 

where p is a constant depending on \vj?\. We remind that the action of SU(2) on the sphere 
S 3 is defined by 

0Ol, Z 2 ) = (0i, 4> 2 )(ZI, Z 2 ) = (4>iZ! - (j) 2 Z 2 , (p!Z 2 + (f) 2 Zl). 

The initial velocity vector v at the point (1,0) has the form v = (i\v e ia ), where we 
write e ia for the initial horizontal velocity v%. According to Remark [3l the coordinate w\ 
can be written as 

Wl = e ~ isM Ux cos (sy/l + M 2 ) + - ^ l sin (sy/l + \v T \ 2 j 

V y/l+\v T \ 2 

where (0f)i = i\vjr\<f)i— e _iQ 2 . Since we are interested only in a fixed value of \wi\ = p > 0, 



we can choose two free parameters in = (0 1; 2 ). We set 0i — P- Then \w 2 \ = a/1 — p 2 , 
and we take arg0 2 = a — | in order to make (0f)i pure imaginary and to simplify the 
calculation of |u>i| 2 . Then 



\ Wl \ 2 = p 2 (cos 2 (sy/l + W\ 2 ) + (K 1 '^ ^ r)2 sin 2 (ay/1 + \v?\ 2 ) 



The equality ^^^p - = 1 gives p 2 = (l + (a/1 + \vjr\ 2 — \vjr\) 2 ) 1 . This proves the 
theorem. □ 

5. Boundary value problem and distance on § 2 ™ -1 

Let us find the distance from the point a to a point p in a fiber, i. e., p G ae luJ , u G (0, tt). 

Theorem 4. The Carnot-Caratheodory distance d c _ c (a,p) from a G § 2ra_1 to the point 
p = ae luJ is 



d c _ c (a,p) = y/u(2ir - uj. 
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Proof. We assume that the geodesic j(s, a, Vy_, \v jr|) is parametrized by arc length and 
u G (0, 7r). We need to solve the equation 

e~ isM (acos (sy/l+ M 2 ) + -^JL==sm ( Sy /l + M 2 )) = ae iu . 
Arguing as in Theorem [2} we conclude that the geodesic intersects the fiber at the moments 



iirn 1 



s = , = n at the points ae v /l+l "^' 2 , n = 1, 2, We are interested in finding 

such that the geodesic 7(5, a, t>%, meets the point ae*^ for the first time, or in 
other words for n — 1. Thus, 

fi l^l A i 12 ( 7r ~ a; ) 2 * V 

7T 1 = = CJ ==>• t>F = — ; r =>• S = W UJ[27r — u). 

\ yiT^p/ 1 • F| Lo(2n-u) V 1 ; 

Since the geodesic is parametrized by arc length, the value s gives us the length of the 
geodesic joining the points a and p. If the point p tends to a (or in other words u — > 0), 
then the velocity \vj\ tends to infinity and the length tends to 0. 

The geodesic is not unique. Varying the directions of the horizontal velocities we obtain 
uncountably many geodesies parametrized by the (2n — 3)-sphere. But all of them have 
the same length. 

There are geodesies that end at point p at times s = , = n, n > 1. The initial 

Vi+M 2 

velocity for the n-th case is 



7m- u 

(M)n 



u(2rm — lu) 



For any fixed u we have 



(\vM)n — > TT - ' $n = \/uj{2ixn - u) — > 

n— >oo ZUl n— >oo 



OO. 



If 00 G (71, 27t), then we can switch u to —u by spherical symmetry. If u G (0, 2tt) + 27rm, 
m — 1, 2, . . ., then in order to find \vt\ we have to solve the equation ( 1 , = ) = UJ+nm 

for different combinations of m and n. The condition UJ ^ m G (0, 1) reduces this case to 
the one considered above with the argument lo G (0, it). 

Since \J u{2ti — a;) is the minimal length among all the geodesies, it gives the Carnot- 
Caratheodory distance, see [20J. □ 



5.1. Boundary value problem and distance on S 3 . Since the study of boundary value 
probem for arbitrary points of § 2n_1 , n > 2 is rather difficult, we concentrate our attention 
on the case of S 3 . As it was shown in [13] the sub-Riemannian structure given by the 
Hopf fibration and by the group structure SU(2) on S 3 coincide. Therefore, we simplify 
considerations taking the base point a = 1 = (1,0) in complex coordinates. If v = 
(vq + iv i,V2 + W3), then vq = 0, v\ = \vj?\, and v%, V3 are arbitrary. The formulas for 



HOPF FIBRATION. 



13 



geodesies turn into (zi(s) , z 2 (s)) , where 



Zl (s) = e~ iVlS (cos(\\v\\s) + i-p- sin(||v||s)), 



z 2 (s) = e 1 — — - — sm( \\v si 



\v | 

In what follows we do not assume parametrization of geodesies by arc length, we suppose 
that all of them are parametrized in the interval [0, 1]. For convenience, we rewrite geodesic 
equations with the following notations 

u := Ti — iT? re := — n — n — > P := M ■ 

IMI \\ v \\ 

Then 

(7) z x (s) = e~ lups (cos ps + iu sin ps) , 

(8) z 2 {s) = re~ i{ups+a) sin ps. 

Let us denote the endpoint at s — 1 of the geodesic by (zi, z 2 ), Z\ = l^ile* 6 * 1 , z 2 = \z 2 \e l62 . 
Then, 

(9) z\ = e~ lup (cos p + iu sin p), 

(10) z 2 = re- i(up+a) sin p. 

Given an endpoint (zi,z 2 ), i.e., given the values of \zi\, arg^i e [— and aigz 2 e 
[— 7T,7r), the unknown parameters are u, p, and a. The last parameter a is the simplest 
one, which one defines it at the end of all computations relating it to axgz 2 . 

Remark 4. Observe that u may belong only to the open interval u G (—1, 1). If ' u = ±1, 
then r = by u 2 + r 2 = 1, and z 2 (s) = 0, which implies that formulas (TpESP reduce to the 
fixed point (1,0). 

• Exceptional cases. Let us start with the cases when the endpoint (zi,z 2 ) lies on 
the vertical line or on the horizontal sphere § 2 . In the first case z 2 = 0. 

— Vertical line and loops. The general case was considered in the previous 
subsection in Theorem |H 

— Antipodal point (—1, 0) is the intersection of the vertical line starting at (1, 0) 
and S> 2 considered as the base space for the fiber at (1, 0). We distinguish this 
case because this point is the intersection of the horizontal sphere S 2 and the 
vertical line, i.e., these points can be connected with geodesies either lying on 
§ 2 for all s e [0, 1], or leaving § 2 . 

Proposition 1. The geodesies connecting the antipodal points (1,0) and (—1,0) 
are given by formulas ([7]) and ([8]), s G [0, 1] with 

* p = Tim, in G N; 

* u = (2p + l)/m for even m, where p is integer and — y < p < y — 1; 

* u = 2p/m for odd m, where p is integer and < p < ^T^i 
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* a is arbitrary. 



The length of these geodesies is d = u^jm 2 — (2p + l) 2 for even m and d = 
7T\/ m? — 4p 2 for odd m. The Carnot-Caratheodory distance is d c _ c = it and it 
is realized by the geodesic with u = 0, p = tt, which lies on S> 2 . 

Proof. The value r ^ 0, see Remark |H So we have p = ttitl > 0. We avoid 

the case m = because the speed is and there is no motion at all. Then 

e~ mmu cos ttitl = — 1 and 

2r> + 1 2r> 

u = — , qeN or u = — -, o6NU{0j, »eZ, 

2q 2q + l 

where m = 2q or 2q + 1 respectively. In the first case, — q < p < q — 1, and 
in the second — q < p < q. The rest of the unknowns are r = \/l — u 2 and an 
arbitrary a. The length of the geodesies is defined as 



= - {2p+ l) 2 , or d 2p , 2q+ i = n^{2q + l) 2 - 4p 2 . 

In the first case, the minimal length of geodesies d = 7r\/3 is realized for q = 1, 
p = or 9 = 1, p = —1. In the second case, d = tt for q = p = 0. Thus, 
the Carnot-Caratheodory distance is given for u = and the corresponding 
geodesic joining (1,0) and (—1,0) lies on S> 2 and it is the half of a big circle 
(mod(a)). □ 

Example. Let us write down explicit formulas for three geometrically different 
(with respect to the rotational symmetry in a) geodesies for u = —1/2, 0, 1/2. 
For u = ±1/2, the geodesies do not lie on S 2 whereas for u = the geodesic 
lies entirely on §> 2 . 

zi(s) = e ±l7TS (cos2irs ± |Msin27rs), J zi(s) ■ 

Z 2 ( S ) = Vl e i(±7r S +Q0 S in27TS, \ z 2 (s) ■■ 

The last geodesic is the minimizer giving the Carnot-Caratheodory distance. 
Horizontal sphere. If the ending point (z\,z 2 ) belongs to the vertical line, 
then there are infinitely many geodesies joining the origin with (£1,2:2). This 
phenomenon is typical for sub-Riemannian geometry and it is seen in many 
examples, for instance, for the Heisenberg group. However, in the case of the 
sub-Riemannian sphere, the number of geodesies joining two points even of the 
horizontal sphere is also infinite. 

Proposition 2. Let (zi,z 2 ) G S 2 . There are countable number of geodesies 
connecting the points (1,0) and (zi,z 2 ). The value Z\ is real. These geodesies 
are given by formulas ([7]) and ([8]), s G [0,1] where 
* P = Pm, itt- ^ N are solutions to the equation 

cos p 



\fzf — cos 2 p 

COS I p — ; ; 

snip 
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/ 2?-COS 2 p m 

* U = U m = ±\ — 

m \l sm z p m 



* a = &rgz 2 -u m p m . 
The Carnot-Caratheodory distance is d c _ c = | arccos zi| and it isrealized by the 
geodesic with u = 0, p = arccos Zi, which lies on § 2 . 

Proof. If the point (21,22) belongs to the sphere § 2 , then Im Z\ = 0. Let us 
assume z\ G (0, 1). This implies 

u sin p cos up = cos p sin up. 

If u = 0, then there is a unique geodesic lying on § 2 joining (1, 0) with (zi, z 2 ) 
modulo repeating big circles, Im z 2 = 0. Its minimal length is d = \ arccoszil- 
Indeed, we have \\v\\ = p in this case which presents the horizontal speed and 
which is preserved under the motion. 

If cos p sin -up = and u 7^ 0, then p = | + rrn or up = irn. In both cases we 
come to the conclusion that u = ±1, which we get rid of, see Remark HJ 
In what follows we consider only the case u G [0, 1) because the case u G (—1,0] 
is treated similarly If cos p sin up does not vanish, then all values sin p, cos p, 
sin up, cos up are non- vanishing too, and hence, up, p 7^ 7rn and up, p 7^ | + Tin. 
So the parameters u and p satisfy the system of equations 
tanup tanp cosp 



up p cos up 

Combining them we obtain the explicit function 

9 zi — cos 2 p 
sin p 

which is defined in each interval p G D n = (arccos Z\ +Tm, 7r(n+l) — arccos Zi), 
which is open because u ^ and we choose arccos z\ G (0,7r/2). Substituting 
u = u(p) in any of two equations from the system (fT2|) we obtain the equation 

cos p 



cos p 



\/ zf — cos 2 p 
I sin p| 



z x . 



Let us denote the left-hand side of the latter equality by <3>(p) for every z\ G 
(—1,1) fixed, see its graph in Figure 1. The function $ is rather complicated 
to investigate completely. 
However, the derivative is calculated as 

$' (p) = _^^ ( l_ u 2 )(1 _ pcotp); 

cos up 

substituting u from (TT5|) . Observe again that the function $ is defined only on 
U^Li D n where it vanishes at the points | + im G D n . We have 

,.7T . (-l) n+1 (l - 2%) . 
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$(p) 




71 




2tt 



3tt 




P 



-1 



-2 



-3 



10 



12 



Figure 1. Graph of the function $(p) for z\ — 0.7. 



The first equation from system (Tl2l) implies that solution to (114j) may be 
searched only for p > tt. The vertical asymptotes are at the points p = p m , 
where p m are the roots of the equation 



P 



cos 2 p 



71 



sin p\ 



+ 7rm, 



m 



1,2,3 



Let us denote by ^(p) the left-hand side of the latter equation. The func- 
tion \I/(p) is continuous in each interval D n and vanishes at its endpoints, see 
Figure 2 and 3. Moreover, 



,7T 



z\{ — h 7vn) < max^(p) < 

2 Dn 



n 



1)tt — arccos^i 



in this interval. 

Thus, there are at least two values of m for n e N, n > n = n (zi), such that 
solutions to equation ffl5|) exist and they are different from | + irn. Let us 
denote by p' n and p" n the solutions to ( |T5l) . such that | + 7rn e (p^, p") and there 
are no other solutions in this interval. Then the function $(p) has asymptotes 
at p' n and p", it is continuous in the interval {p' n ,p'n), vanishes at | + 7rn, and 
ranges in (—00,00) on (p' n ,p'n)- Therefore, the equation (TTJJ has at least one 
solution in (p' n ,p"). Since n ranges in N, n > n for n sufficiently large, we 
have an infinite number of geodesies joining two points on the horizontal sphere 

«2 



The inequality — x 2 > arccosx for x G (0,1) 
entirely lying on S 2 realizes the minimal distance. 



assures that the geodesic 

□ 
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Figure 2. Graph of the function \^(p) for z x = 0.7. 



2i (f + 4vr) 



arccos zi + 47r 




13.5 



14.0 



14.5 



p 



Figure 3. Zoom of the graph of the function \I/(p) for z\ = 0.7 in the interval D4. 



General case. Assume that the endpoint (z\, z%) does belong neither to the vertical 
line nor to the horizontal sphere S 2 . The equation ( flOj) implies 



rp > r\ sin p\ = \z 2 \ 
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-0.4 
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0.4 
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-0.2 














-0.4 











Figure 4. Graph of the function B{u). 
where rp is the length of the geodesic. At the same time, the relation 

COt p = 0\0 2 - 





2 


— U 2 




22 





a x a 2 - 



/ r 2 — I z 2 1 2 



R2 



follows from ([10]) with 

0"i = sgn [cos p] , 0"2 = sgn [sin p] . 

Observe that p = irm corresponds to the exceptional case (z±, z 2 ) in the vertical line, 
and hence, the function cot p is well-defined and finite. The equation implies 



cos p + iu sin p °"i 's/l^i I 2 ~ u 2 + io 2 u\z 2 \ 



\/ cos 2 p + v? sin 2 p 



r\Zi\ 



, p 7rm. 



Case o"! > 0, o"2 > 0. The equations for w and p > become 

27rg, geNU{0}, 



6 1 ! = arccot 



p = arccot 



V\zi\ 2 -u 2 
u\z 2 \ 





12 


— M 2 




22 





m arccot 





2 


— U 2 




2 2 





+ 27r(p — uq), pel*. 



where the branch of arccot is chosen to be in the interval (0, 7r). Observe that 

IT 

arccot ; — ; G (0, —J. 





2 


— M 2 
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;4^-0.2 




-~^~-^Q,2 


0.4 "~~0^^_ 






-10 










-20 







Figure 5. Graph of the function B{u) + 2iv{p — 3u), \zi\ — 0.7. 



If u > 0, then also 



Denote 



U 2 , 71", 



arccot — j — j G (0, 

u\z 2 \ 2 







2 _ 


- U 2 


u\z 2 





B(u) = arccot — : — : u arccot 





2 


— U 2 




22 





One finds its graph in Figure 4. Obviously it is an odd function and |-B(u)| < |. 
The inequalities — n < Q\ < n and \u\ < 1 imply — ^ < 2ir(p — uq) < 
Therefore, — | — q < p < q + ^, or— q < p < q. The equation 9 1 = B{u) 
corresponds to the choice q = 0, p = 0. We visualize possible choices of p for 
q = 3 in Figure 5. 

Case oi < 0, cr 2 < 0. The equations for it and p > turn into 





12 


— U 2 




Z-2 





p = arccot I , h7r(2g+l), geNU{0}, 





Z\ 


2 _ 


- U 2 


u\z 2 







2 


— U 2 
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9\ = arccot — . — . u arccot — — . — . h vr((2p + 1) — u(2q + 1)) 



p G Z with the same branch of arccot as in the previous case. The inequalities 
— 7r < Q\ < 7i and \u\ < 1 imply — q — 1 < p < q + 1. 
— Cases aicr 2 = — 1. These two cases come down to the previous ones changing 
p->-p. 

Carnot-Caratheodory distance. The length of a geodesic is calculated as 



pyl — u 2 , where u and p are defined below. 
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Figure 6. Graph of the function B(u) + 2ir(p - u) for p = -1, 0, 1. 

- Let < \9i\ < — \zx\). In the case <j\ > 0, cr 2 > 0, there is a unique solution 
it = uq to the equation B{u) = 8\ (see Figure 4). The corresponding value for p 
with q = we denote by p . All solutions to the equation B(u)+27i(p— uq) = 9i 
with q > we denote by u VA and the corresponding value of p we denote by 
p Pt q. The function 





12 


— X 2 




22 





Vl — x 2 |^arccot 
is even and increases in x G [0, \zx\], which implies that 

i 



Po\/l — ul = y 1 — Uq I arccot 





12 


-ul 




Z2 





\ Z A 2 - U l,q 

- \l l ~ u pa ( arccot ~ r~j + 2?rg ] = p p , q \/l - u 2 M . 

Similarly other choices of <j\ and o"2 do not give the minimizing geodesic. So if 
< | Q\ | < f ( 1 — | Z\ | ) , then the minimal length is 





12 






^2 





d = y 1 — Uq ^arccot 

where «o is a unique solution to the equation B(u) = 6i. 
Let |^i | > |(1 — |). Then there is no solution to the equation B(u) = 9 1 . 
The situation is even more complicated. We visualize it in Figure 6. If q — 1 we 
can guarantee a solution to the equation B(u) + 2n(p — u) — d\ for \z^\ > 3/4. 
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Therefore, in the case when the equation B{u) + 2ix(p — u) = 9 1 has a solution, 
let us consider q — 1 and q > 1. The function 



Vl — x 2 arccot — : : h 2% 





Zl 






Z2\ 



attains its maximum in the interval x G [— \z\\, |zi|] at the point x — 0. More- 
over, we have an elementary trigonometric inequality 



arccot < 2nx, x £ (0, 1 . 

x 

Hence, in the case o\ > 0, cr 2 > and q > 2, we have the following chain of 
inequalities 



z i\ 2 -<i 



p p ,i\Jl — Up t i = \Jl — it 2 ! J arccot j — j — — 1 < arccot y^-r < 



zA 2 - ul 



< 2nq\z 2 \ < yj 1 - u 2 pq I arccot ^ — — + 2vrg = Pp ^l-u 2 pq . 



Thus, d = min p v \\l\ — ul-,. 

p=-i,o,i V p ' 
In general case, we are able to give an algorithm of finding geodesies and 

the length. Let q m be a minimal non-negative integer for which the equation 

B(u) + 2n(p — qu) = 9\ has a solution. Then the distance is calculated as 



rf = min|p Pi9m y / l- M 2^ 



among p = {-q m , -q m + 1, . . . , -1, 0, 1, . . . , q m - 1, q„ 
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